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Abstract. The many-body Green's function theory developed in our previous 
work for treating the reorientation of the magnetization of thin ferromagnetic films 
is extended to include the exchange anisotropy. This leads to additional momentum 
dependencies which require some non-trivial changes in the formalism. The theory is 
developed for arbitrary spin values S and for multilayers. The effects of the exchange 
anisotropy and the single-ion anisotropy, which was treated in our earlier work, on 
the magnetic properties of thin ferromagnetic films are compared. 
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1. Introduction 

There is increasing activity in experimental and theoretical investigations of thin 
magnetic films and multilayers. Of particular interest is the reorientation of the 
magnetization as function of temperature and film thickness. The simplest ap- 
proach for treating thin ferromagnetic films is the application of mean field theory 
(MFT) to a Heisenberg model, by either diagonalizing the corresponding single- 
particle Hamiltonian |lj or applying thermodynamic perturbation theory (2j. This 
approximation completely neglects collective excitations (magnons = spin waves). 

In order to take the influence of these collective excitations into account, we have 
turned to many-body Green's function theory (GFT), which allows reliable calcu- 
lations over the entire range of temperature of interest. In reference P] we treated 
a spin S = 1/2 Heisenberg monolayer in a magnetic field and found, by compar- 
ing with Quantum Monte Carlo (QMC) calculations, that a Tyablikov (RPA) |3] 
decoupling is a very good approximation. Therefore, we did not try to go beyond 
RPA in the subsequent paper [5], in which we treated the field-induced magnetic 
reorientation of a ferromagnetic S = 1 monolayer, whereby a second-order single-ion 
anisotropy was also included. Whereas the exchange interaction terms are decou- 
pled by RPA, this is not allowed for the terms coming from the single-ion anisotropy 
because this leads to unphysical results. Instead, we applied a decoupling procedure 
proposed by Anderson and Callen [U], which, however, is a good approximation only 
for small anisotropics. This was shown to be the case in reference |f7], where we 
were able to treat the single-ion anisotropy exactly (for any strength) by introduc- 
ing higher-order Green's functions and subsequently taking advantage of relations 
between products of spin operators, which leads to an automatic closure of the hi- 
erarchy of the equations of motion for the Green's functions with respect to the 
anisotropy terms; the terms from the exchange interaction are still decoupled by a 
generalized RPA scheme. In reference PJ we have investigated the quality of this 
approach by comparing with QMC calculations. In reference 17] we treated the spin 
5 = 1 case only; the formal generalization to spins 5* > 1 is possible, but its nu- 
merical realisation is quite cumbersome. This is not the case when remaining at 
the level of the lowest-order Green's functions and applying the Anderson- Callen 
decoupling. In this case, not only is the treatment of spins 5* > 1 feasible, but also 
multilayers can be described, as was done in reference P|. To make the treatment 
of multilayers practicable, we had to apply a new method which not only uses the 
eigenvalues but also the eigenvectors of the non-symmetric matrix which governs 
the equations of motion for the Green's functions. We mention a few related pa- 
pers in which Green's function theory is also applied to spin reorientation problems. 



2 



In reference jTOj Green's functions are applied to the reorientation problem after 
a Holstein-Primakoff mapping to bosons, which is only a valid description at low 
temperatures. In reference |TJ the reorientation is obtained by the competition of 
an approximately decoupled single-ion anisotropy and a schematic shape anisotropy. 
In the interplay of the dipole coupling with an easy-plane single-ion anisotropy 
is discussed, and in [13^ the in-plane dipole couphng anisotropy of a square ferro- 
magnetic Heisenberg monolayer is considered. 

In our previous work mentioned above, we treated an isotropic Heisenberg ex- 
change interaction plus a single-ion anisotropy, the magnetic dipole coupling, and an 
external magnetic field. In the present paper, we include the exchange anisotropy in 
the formalism and discuss its effect on the magnetic properties of thin ferromagnetic 
films, in particular in comparison to the influence of the single-ion anisotropy. As 
will be shown, the treatment of the exchange anisotropy requires some non-trivial 
changes in the formalism, which are necessary due to additional momentum depen- 
dencies which are absent in the treatment of the single-ion anisotropies. 

The paper is organized as follows. In Section 2 we establish the Green's function 
formalism. For pedagogical reasons, a large part of the formalism is demonstrated 
for the monolayer case and subsequently generalised to the multilayer case, which 
can easily be done. In Section 3, we discuss the field-induced reorientation of the 
magnetization when using the exchange anisotropy for determining the orientation 
at temperature T = 0. We fix the strength of the exchange anisotropy in such a 
way that it produces the same Curie temperature for a monolayer as the application 
of the single-ion anisotropy strength used in our previous work E] . This enables 
a comparison between the roles of the single-ion and exchange anisotropies on the 
magnetic properties of thin ferromagnetic films. Section 4 contains a summary of 
the results. 

2. The Green's function formalism 

We formulate the theory in such a way that the results of our previous work [SJ E] 
are obtained as limiting cases. 

We consider a spin Hamiltonian consisting of an isotropic Heisenberg exchange 
interaction with strength Jm between nearest neighbour lattice sites, an exchange 
anisotropy with strength Dki, a second-order single-ion lattice anisotropy with strength 
K2^k, a magnetic dipole coupling with strength g^i, and an external magnetic field 
B = {B'',By,B'): 

n= - W Ms^s^ + s^sn - ^ E DMS'kS[ - E K2AS'k? 

^ <kl> ^ <kl> k 
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^ kl ^kl 



(1) 



Here the notation = ± iS^ and = B^ ± iB^ is introduced, where k and 
I are lattice site indices and < kl > indicates summation over nearest neighbours 
only. The only difference from reference jH] is the additional exchange-anisotropy 
term. 

In order to treat the reorientation problem for general spin S, we need the 
following Green's functions 



Ga,mn/ \ / / oa. f Qz\mf Q—\n\\ 



(2) 



where a = (+, — , z) takes care of all directions in space, rj = ±1 refers to the 
anticommutator or commutator Green's functions, respectively, and n > l,m > 
are positive integers, necessary for dealing with higher spin values S. 
The exact equations of motion are 



Gci,mn / \ 



with the inhomogeneities 
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(3) 



(4) 



where (...) = Tr{...e ^^). The equations, leaving out for the moment the terms due 
to the dipole couphng, are given explicitly by 
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(5) 



After solving these equations the components of the magnetization can be de- 
termined from the Green's functions via the spectral theorem. A solution is pos- 
sible by establishing a closed system of equations by decoupling the higher-order 
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Green's functions on the right hand sides. Contrary to reference j7], where we pro- 
ceed to higher-order Green's functions, we stay here at the level of the lowest-order 
equations. For the exchange-interaction and exchange-anisotropy terms, we use a 
generalized Tyablikov- (or RPA-) decoupling 



•/3,mn 



(6) 



The terms from the single-ion anisotropy have to be decoupled differently, because 
an RPA decoupling leads to unphysical results; e.g. for spin 5* = 1/2, the terms 
due to the single-ion anisotropy do not vanish in RPA, as they should do, because 
in this case Si -^2,j(('S'f )^) is a constant and should not influence the equations of 
motion. In the appendix of ref. [H] we investigated different decoupling schemes 
proposed in the literature, e.g. those of Lines [T3| or that of Anderson and Callen 
[H], which should be reasonable for single-ion anisotropies small compared to the 
exchange interaction. We found the Anderson- Callen decoupling to be most ade- 
quate. It consists in implementing the suggestion of Callen ^3] to improve the RPA 
by treating the diagonal terms arising from the single-ion anisotropy as well. This 
leads to 



{{isrs!+s!sry,isxisir))r 



2S- 



:[S{S+l)-{S!S:)]]Gp-. 



(7) 



This term vanishes for S* = 1/2 as it should. 

After a Fourier transform to momentum space, one obtains, for a ferromagnetic 
film with N layers, 3A^ equations of motion for a 3A^-dimensional Green's function 
vector G"^^: 



(8) 



where 1 is the 3A^x3A^ unit matrix. The Green's function vectors and inhomogeneity 
vectors each consist of N three-dimensional subvectors which are characterized by 
the layer indices i and j 



G-(k,a;) 



Gtri^,uj) ) 
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mn 
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(9) 



1 /,z,mn 1 



The equations of motion are then expressed in terms of these layer vectors, and 
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3x3 submatrices Tij of the 3A^ x 3A^ matrix F 



col 
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j = l,...,N. 



(10) 

After applying the decouphng procedures (jH)) and ((Tj), the F matrix reduces to a 
band matrix with zeros in the Fjj sub- matrices, when j > i + 1 and j < i — 1. The 
diagonal sub-matrices Tu are of size 3x3 and have the form 



Hi 








Hi 




'11^ 



where 



Ht 



Z, + {St).Mq-j^) , 

+ Duq{S^) + {Ju+, + A,i+i)(5'^i) + iJ^,^-l + A,-l)(^tl) 

+K,,2{Sn{l-^,[S{S+l)~{S:sn]) , 

Bt + {St)Mq - 7k) + J^,^+l{St+,) + J^,^-l{Stl) , 



121 



For a square lattice and a lattice constant taken to be unity, 7k = 2(cos kx + cos ky), 
and g = 4 is the number of intra-layer nearest neighbours. Except for the exchange- 
anisotropy terms, Dij, these equations are the same as in reference |9:. Putting all 
Dij = 0, one has Hf = Hf. Note that owing to the momentum dependence in 
coming from the exchange anisotropy, Hf ^ Hf, which forbids a naive extension 
of the formalism of reference j^, as discussed below. 

Approximating the dipole coupling by mean field theory (MFT), which is a 
good approximation when the dipole coupling strength is small as compared to the 
exchange interaction strength ( as proved in appendix A of 0), one finds that the 
effects of the dipole coupling can be included as an effective field: 



TV 



B~- 



± 



N 



Bl = B^-2j2g,,{S;)T 



where the lattice sums for a two-dimensional square lattice are given by (n 
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n 



Im 



(/2 + m2 + ^2)5/2- 



(13) 

(14) 
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The indices Im run over all sites of the j*'* layer excluding terms with P+m^+n'^ = 0. 

One observes that the dipole coupling in MFT leads to a renormalization of the 
external field: there is an enhancement of the transverse fields, and a reduction of 
the field perpendicular to the film. 

The 3x3 off-diagonal sub-matrices Fij for j = i ± 1 are of the form 







-(j,, + A,)(5r) 





(15) 



J,,(5f) 

When treating the monolayer, one can use the spectral theorem for calculating 
the components of the magnetization. This was done in reference j5] for the case of 
spin S = 1 and the single-ion anisotropy by using the commutator Green's functions. 
In order to obtain sufficient equations it was necessary, to add equations coming from 
the condition that the commutator Green's functions have to be regular at = 0, 
which we call the regularity conditions. 

For the multilayer problem, however, a naive application of the spectral theorem 
turned out to be forbiddingly difficult. Instead we invented a method, which we call 
the eigenvector method, that uses the eigenvectors as well as the eigenvalues of the 
F-matrix governing the equations of motion. This opened up a practicable way to 
treat multilayers [H]. 

If instead, anticommutator Green's functions are used, it is not necessary to 
introduce the regularity conditions, which nevertheless are valid. We demonstrate 
this explicitly for the monolayer. The use of the anticommutator Green's functions 
also suggests a way for finding the procedure which can deal with the additional 
k-dependencies coming from the exchange anisotropy. 

For simplicity, we consider the reorientation in the x — 2;-plane, i.e. we use as 
external field B = (5^,0,5^). The equations of motion for the monolayer in this 
case are 



/ 
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(16) 



This system of equations has three eigenvalues 



(17) 



and the equations are solved by 



\a,mn 
V 
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where A"'™-" is the determinant of the matrix in equation where column a is 
replaced by the inhomogeneity vector, and A = uj{uj — e\^{uj + Ck). 
Now the spectral theorem ^H] is applied in momentum space 

= {{S^T{S~rS-\ = lim^ f ^^(G™ 

(19) 

Using the relation between anticommutator and commutator 

where it is important that the commutator inhomogeneities A'^^i do not depend 
on the momentum k, one obtains the following set of equations 

^.^+,mn _ ^.^^,^n ^ ^Jrnn ^ 1 ^^th(^) - ^H^) + A^^^^, (21) 

_^.^-,mn ^ jj^c^n ^ ^-^«( COth(^) + ^H^ - (22) 
^,,^+,mn _ ffX(j-,mn ^ ^ffx ^^-^_n _ ^+^_n^ _ coth(^) . (23) 



Because the observable correlations are in real space, we have to perform a corre- 

1 ^a,mn 



sponding Fourier transformation {{S^)^{S )"5'f) = Cf'™" = jrJ2kC] 



Fourier transform of equation (|22p yields 

_c-'"^"+- y f^c''"^^ = -A-i"'i+-A-f'i- y ^ coth(^)--A^:!""i- y — . 

(24) 

Putting this into the Fourier transform of equation ()21|). one can eliminate the term 
jf Y^k ^C^'™"". This will turn out to be important in the later discussion of the 
eigenvector method where the formalism has to be modified because one cannot 
take the k-dependent terms outside the sum (integral). One obtains 

C^'^-Cp^-^-iA-r^-A^^--) = ^(A-i!!1 + <=f^i)^E;^coth(^). (25) 

The Fourier transform of equation ()23|1 can be done directly and gives 

- Cp- - \{A-fP[ - At,fP[) = E coth(^). (26) 

Equations and are sufficient to determine the observables. 

To elucidate these equations we derive the exphcit expressions for spins S" = 1/2 
and S = 1. 

For 5* = 1/2 we need m = and n = 1. This gives two equations of motion 
which determine (S^) and (S^). 
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From equation one finds with {S^ S^) = 1/2 — {Sf) and {S^ S,^ ) = and 



A 



,01 



- and A;:!'_, = 2{S!) 

1 



- = (^-)lyiiLcoth(^^ 



and from equation witli A^'_i = —{S^ ) = — (5'f) 



1 = (5^)lyiLcotli(^^ 

2 ^ ^ N ^ H '^ ^ 2 ' 



(27) 



(2^ 



Tliis are two equations wliicli determine tlie two unknowns (S^) and (S^). No 
regularity conditions are necessary. 

For 5 = 1 one needs equations (j25p and (j26p not only for (n = 1, m = 0) but also 
for (n = l,m = 1), (n = 2,m = 0), (n = 3,m = 0). This yields 8 equations for the 
eight unknowns {S~), (S^), {S-S'), {S'S~), {S'S'), (S'S'S-), (S'S'S-), {S-S-S'S'). 

The left sides of equations (|25p and (|26p are the same. Therefore we write them 
pairwise. For (n = 1, m = 0) we have 



(5^)^Ek|lcoth(f^) 



2-{S'S') - {S-S-) 



(29) 



(5-)iEk|7Coth(f^). 
For (n = 1, m = 1) and {S^S^'S^} = (5*^), valid for 5* = 1, we have 

h{S' + {S'S') - {S-S-) - 2{S'S-S-) - 2) = 
' \[{S-S-) + 3(5^5-) - (5^) - 2) ^ Ek coth(f^) 

< 

. (5^5-)iEk|rCoth(%^). 
For (n = 2, m = 0) and (S^S^S^) = we find 

3{S-) - {S-S') - {S-S'S') + 2{S'S~) = 
(2{S^S-) + {S-))j^E^j^coth{Sf) 

2(^-^-)^Ek|7Coth(f^). 
For (ra = 3, m = 0) we have with {S^ S^) = and (S^S^S^) = 

2(5-5-) - {3-3- S') - {3-3- 3'3') = 
' (3(5^5-5-) + 3(5-5-))iEkf|coth(%^) 

0. 



(30) 



(31) 



(32) 
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From equations (021) we find {S-S-S'S') = 2{S-S-) - {S-S-S') and {S'S-S-) = 
— {S^S^). The remaining correlations are determined by the previous six equations. 

Instead of using these 8 equations one can also apply the regularity conditions, 
as we did in reference jS] when working with the commutator Green's functions. 
Then one can express all correlations in terms of the correlations (S^) and (S^S^) 
in the case of spin S = 1, and in terms of the moments ((5*^)") with (n=l,..., 2S) 
for arbitrary spin S, and one need only solve two equations in the S = 1 case or 25 
equations for arbitrary spin values 5*. 

The regularity conditions are obtained from the fact that the commutator Green's 
function has to be regular at the origin 



lima;Gk,T=-i = 0' (33) 



which leads to the relations 

H^'A^^I + H'^A-fT^ + 2iyM;f!"i = 0. (34) 

Note that these relations are also obtained by equating equations and (j^HI)- 
For m = 0, n = 1 we obtain the first regularity condition 

{S^) = ^{S^y, (35) 

i.e. the knowledge of (S^) determines (5"^). 

From the definitions (fT^ one sees that the prefactor is momentum independent 
and the relation generalized to the multilayer can be written as 

Hf Zi 

In the case of spin 5* = 1, equation (j25|l with [n = l,m = 0) and {n = l,m = 1) 
together with the regularity conditions determines all desired correlations for spin 
S = 1. We demonstrate this by deriving the equations for {S^) and {S^S^) already 
derived with the commutator Green's functions in reference jS]. 

Equation ^ gives for (n = 1, m = 0) with {S-S+) = 2- {S') - {S'S') 



1 / WW Be^ 
2 - {S^S^) - {S-S-) = (5^)- E V 1 + 777777 coth(^). (37) 



iV V V WH' ' 2 

k 

Prom the regularity conditions ()34|) one finds for S=l 

(S-S-) = — (3{S'S^) - 2). (3^ 

Putting this into equation (jH7|) gives the first of the desired equations 



4 - 2{S'S')(1 + i^r) - iS') (2 - (|;)'^) ^ E V IH- coth(^|i) = 0. (39) 
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The second equation is obtained from equation (j25j) with [n = l,m = 1), and the 
regularity conditions which relates (S^S^S^) to {S^) and {S^S^). This leads to 



{S'){2 - (f!)^) - 2{3{S'S') - 2)15:^1 + |!§;coth(§i) = 0. (40) 

k 

The only difference from the corresponding equations of reference is that the 
square root cannot be taken out of the sum (integral over the first Brillouin zone) 
because of the momentum dependence of coming from the exchange anisotropy. 
If the latter is zero as in reference [S], only the single-ion anisotropy survives and 
one recovers the original equations. 

The explicit derivations above are done for pedagogical reasons, but they would 
not have been necessary because one obtains by subtracting equation (PT|l from 
equation 



liA-fll - coth(^) + + A-r-1)el§ 



X 



which corresponds to equation (27) of reference [S] which is the starting point for 
deriving the equations for the moments explicitly. 

For the treatment of multilayers we have to use the eigenvector method as men- 
tioned above. The essential features are as follows. One starts with a transformation, 
which diagonalizes the F-matrix of equation (jH} 

LTR = n, (42) 

where f2 is a diagonal matrix with eigenvalues Ur {t = 1, 3N), and the transfor- 
mation matrix R and its inverse R^^ = L are obtained from the right eigenvectors 
of r as columns and from the left eigenvectors as rows, respectively. These matrices 
are normalized to unity: RL=LR=1. 

Multiplying the equation of motion (jH)) from the left by L and inserting 1=RL 
one finds 

(cui - o)lg;;^" = la;;^'^. (43) 

Defining = LG;;^" and = LA™" one obtains 

{ul - ri)e?7" = ^7". (44) 

Ql^"" is a new vector of Green's functions, each component r of which has but a 
single pole 

Am,n,T 

g^n,r ^ _Jl _ ^45) 
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This is the important point and allows application of the spectral theorem to each 
component separately. This gives with C™" = LC^ 



^mn,r ^ + ^(1 - T])- lim ''=+\ (46) 

In reference we used the commutator [t] = —1). Here, we proceed with the 
anticommutator (r^ = +1), so that the second term in equation ()4fjj) is zero and one 
obtains the original correlation vector C™" by multiplying C™" from the left with 
R; i.e. 

C™" = RSLA^^^^, (47) 

where £^ is a diagonal matrix with matrix elements Sij = 5ij{e^^^ + 1)~^. With the 
relation we find 

f^mn ^ K£L(Ap_i + 20'"'^) , (48) 

or 

C™" = (1 - 2Rf L)-iRf LAp_i. (49) 

For the monolayer it can be shown explicitly that the eigenvector method yields 
the equations ()2ip22|23j) derived before. In this case the eigenvectors by which the 
transformation matrices R and L are constructed can be given explicitly. They are 



R 



V 1 1 1 J 



(50) 



and 



4eE 



-{ek + H^)H'' {ek-H')H'' 2i7^^^ . (51) 

Putting this into equation yields equations (j21l22l23|) . 

In order to obtain the correlations in real space, equation (jl^ has to be Fourier 
transformed and the resulting integral equation has to be solved self-consistently. 
By inspecting the expressions for the monolayer, one can show that the inverse 
(1 — 2R£^L)~^ does not exist . Therefore this equation cannot furnish the solution. 
However, one can show that the formulation with the anticommutator relation can 
be transformed into the result for the commutator relation: £^ is a diagonal matrix 
with Sij = 5ij{e^^^ + l)-i. With the relation (R^L)-i = L-i(^)"^R"^ = R^'^L 
one obtains from equation (j49p 

C = (R(l - 2f )L)"iR^LA 
= R(1-2^)"^^LA 

= R£LA, (52) 
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where 8ij = Sij ^^^e ~ ^iji^^"^' ^' which still is of no use because it diverges for 
Ui = 0. 

In reference jH], it was shown that the matrix RqLo, where the index refers to the 
eigenvectors with eigenvalue zero, is a projection operator onto the zero eigenvalue 
space. The situation can then be remedied by projecting equation (f^ onto the 
no n- zero eigenvalue space with 1 — RqLo, which leads to the commutator expression 
for the correlations of ref. [H] 

C = R£oLA + RoLoC, (53) 

where Sq is equal to S in which the diagonal elements corresponding to u; = have 
been set to zero. 

The problem one is confronted with in applying the eigenvector method to this 
equation is that the the exchange anisotropy introduces a momentum dependence 
in the projection operator. Then the projector cannot be taken out of the integral 
when a Fourier transform to real space is performed as in the case of the Anderson 
Callen decoupling of the single-ion anisotropy only. The way out is to eliminate the 
projector by a transformation in one component of equation ()53|) . which is sufficient 
to establish the integral equations of the eigenvector method. This procedure was 
inspired by the elimination of the disturbing term in equations (j21|22j) . 

The adequate transformation is found to be 



/ 1 


1 


0^ 




(l 


-1 o\ 


-1 


1 





T = 


1 


1 


V 





2^ 




lo 


l) 



with T-^T = 1. 

Applying this transformation to equation 

T-^C = T-^R^oLTT"^A + T-^RoLoTT-^C (55) 

and inserting the eigenvectors ()50I51|) transforms the second component of the vector 
T-iRoLoTT-^C to zero, and one recovers equation from the second row of the 
transformed equation (|55p . from which, together with the regularity conditions, the 
integral equations for the correlations for each (m, n)-pair are obtained. 

The eigenvector method is then immediately generalized to the case of N layers 
by applying the transformation to equation ()53|) read as a 3A-dimensional problem, 
thus constructing 3N x 3A-matrices with sub-matrices formed with the transforma- 
tion (j54j) on the diagonals. 
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3. Numerical results 



In this section we discuss results of calculations for a square lattice including the 
exchange anisotropy. The integral equations ()55p for determining the correlations 
together with the regularity conditions as derived in Appendix A of this paper are 
solved self-consistently by the curve following method, which we described in detail 
in Appendix A of reference |T|. A comparison with the results for the single-ion 
anisotropy as used in our previous work [SJ |H] is effected by fixing the strength of 
the exchange anisotropy such that the Curie temperature of a square monolayer 
with spin S = 1 agrees approximately with that of a single-ion anisotropy with 
strength K2 = 1, which was used in most of the calculations of our previous work. 
The exchange interaction strength, J = 100, and the dipole coupling strength, 
g = 0.018 (corresponding to the case of Ni), are taken to be the same in both kinds 
of calculations. The exchange anisotropy coupling strength turns out to be D=0.7. 
In order to compare results for different spin values all parameters are scaled as 
J J/S{S + l),g^ g/S{S + 1),D ^ D/S{S + 1), K2 K2/S{S - 1/2) and if 
a magnetic field is applied, B — > B/S*. 



Magnetization in z-direction: S=1 J=1 00 
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Figure 1: The magnetization {S^) and {S^S^) of a ferromagnetic spin S = 1 Heisen- 
berg monolayer for a square lattice are shown as functions of the temperature (no 
magnetic field). Comparison is made between Green's function (RPA) calculations 
using the exchange anisotropy {D = 0.7, open circles) and the single-ion anisotropy 
{K2 = 1, solid dots) with Anderson-Callen decoupling. The corresponding results 
of mean field (MFT) calculations are also shown. 



In figure 1 we display the magnetization (S^) and its second moment (S^S^) 
as functions of the temperature for a S* = 1 monolayer using Green's function 
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theory (denoted as RPA). Using a value of D — 0.7 as the strength of the exchange 
anisotropy yields nearly the same curve as for an Anderson-Callen decoupling of the 
single-ion anisotropy with a strength of K2 — 1. The results for the corresponding 
mean field (MFT) calculations also practically coincide with each other. As is well 
known and also discussed in our previous work, the neglect of magnons results in a 
Curie temperature which is more than a factor of two larger than that obtained by 
including magnon excitations. This difference for a monolayer is much larger than 
the corresponding difference for bulk ferromagnets. 



reorientation in magnetic field at fixed temperature 
S=1 J=1 00 g=0.066 D=0.7 T=30 
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Figure 2: The components of the magnetization (S^) and (S^) and the absolute 
value iS at a fixed temperature T = 30 as function of an external magnetic field in 
the x-direction, B^, are shown for a ferromagnetic spin S = 1 Heisenberg monolayer 
for a square lattice. Also shown are the equilibrium reorientation angle and the 
critical reorientation field, Bj^, at which in-plane orientation is reached. 



In figure 2 we show the reorientation of the magnetization at a fixed temperature 
due to a transverse field in x-direction, B^. The magnetic field in z-direction is set 
to zero, B^ — 0. In this case, the strength of the dipole coupling is chosen to be 
g — 0.066 (a value corresponding to Co). As a function of the external field, the 
x-component of the magnetization < 5"^ > rises linearly, whereas its 2;-component 
< iS^ > falls to zero, where the in-plane magnetization is reached (^^o = 90°). The 
absolute value S = ^/ < >^ + < >^ remains constant, as it should do, and is 
also shown in the figure. 

In figure 3, we show the magnetization components as functions of the temper- 
ature for the same parameters as in figure 2 at a constant external field, B^ = 0.3. 
The component (S^) stays constant until the component (S^) has dropped to zero, 
and an in-plane magnetization (6*0 = 90°) is reached. The fact that {S^) is constant 
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reorientation as function of the temperature 
at fixed B''=0.3 J=1 00 g=0.066 D=0.7 S=1 




temperature 

Figure 3: The components of the magnetization (5*^) and (5*^) and its absolute value 
S for a fixed magnetic field = 0.3 as function of the temperature are shown for 
a ferromagnetic spin S = 1 Heisenberg monolayer for a square lattice. Also shown 
are the equilibrium reorientation angle 6q and the critical reorientation temperature, 
Tr, at which in-plane orientation is reached. The small horizontal arrow indicates 
the value of (5*^) below which complex eigenvalues occur. 



for temperatures below the reorientation temperature, T < Tr, can be understood 
from the regularity condition ()35|1 . which, for a monolayer with scaled parameters 
{{S^) drops out of the expression if = 0), is (S^) / S = (^4£)_^^oys{s+i) • "^^^^ 
expression restricts the range of parameters for the validity of our approach be- 
cause {S^)/S must be < 1. Above the reorientation temperature, (S^) remains zero 
whereas (S^) remains finite because of the field in the x-direction and decreases 
slowly with increasing temperature. The absolute value of the magnetization S is 
also shown; above Tr one has S = (S^). 

With the chosen parameters one observes a novel feature in the calculations 
which is connected with the introduction of the exchange anisotropy together with 
the dipole coupling. In this case, the eigenvalues of equation (fTTj) become complex 
above a certain temperature, i.e. below a certain value of (5*^). This behaviour 
occurs quite naturally in the theory. Because the F-matrix of equation (jH)) is real, 
its eigenvalues and eigenvectors, if complex, occur pairwise as complex conjugates, 
and the term R£^L in the equation (j55p is real, so that one has to do with a real 
integral equation. The complex eigenvalues and vectors have to be taken seriously, 
and are necessary for obtaining the results of figure 3. The complex eigenvalues are 
connected with the additional term in Hf of equation (fT^ . which comes from the 
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exchange anisotropy. This can be seen analytically by considering the dispersion 
relation ()17|). which shows that the optimal condition for the occurrence of complex 
eigenvalues is at k = ( thus 7k = 4). With q=4, = 0, K2 = one has for the 
monolayer from equation (jl7|) 



-k=0 



(56) 



Complex eigenvalues occur if the second term is larger than the first term. Using the 
regularity condition (jHTj) for (S*^), one obtains (using scaled parameters for 5 = 1) 
a condition for the occurrence of complex eigenvalues, 



25^ 



2gT^ 



{AD - ?,gT^) \ {AD - 2gT^) ' 



(57) 



This yields complex eigenvalues for (5*^) < 0.1157 at T ~ 53 for the parameters 
used in figure 3. This is confirmed in the numerical calculations, where of course 
complex eigenvectors also occur at finite k. Enlarging B^ and/or g increases the 
range of complex eigenvalues. No complex eigenvalues occur if the dipole coupling is 
set equal to zero. In our previous work where we used the single-ion anisotropy 
together with the dipole coupling, complex eigenvalues never occured, which can 
be understood from the structure of the dispersion relation (fTTj) by putting the 
exchange anisotropy terms equal to zero. 
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Figure 4: Normalized magnetization curves for {S^)/S and {S^)/S for a monolayer 
with spin values between S = 1/2 and 5* = 6 shown as functions of the temperature 
using the exchange anisotropy. The corresponding results obtained with a single-ion 
anisotropy {K2 = 1) are shown in the inset (from reference j[9]). 
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In Fig. 4 we display the normalized magnetizations {S^)/S of a monolayer as 
functions of the temperature T for all half-integral and integral spin values ranging 
from 5' = l/2to5' = 6. In the case of the exchange anisotropy, the reorientation 
temperature Tr is practically the same for all spin values. One observes a spin 
dependence of the magnetization curves which decreases with increasing spin. The 
curves, however, saturate very quickly; on the scale of the figure, one cannot distin- 
guish the curves between S — 2 and S — 6. When using the single-ion anisotropy 
(the results are shown as an inset) there is a difference in the reorientation temper- 
atures, but the magnetization curves again saturate very quickly (approaching the 
classical hmit) but in the opposite direction. In both cases, the values for {S'-^)/S 
remain rather small owing to the application of the small field in the x-direction 

= 0.1. 




Figure 5: Sublayer magnetization components (Sf) and (Sf) as functions of the tem- 
perature for spin S — 1 films with N layers calculated with the exchange anisotropy. 
The reorientation temperatures can be read from the curve in the N-T plane. 

In figure 5, we display the sublayer magnetization components (Sf) and (Sf) 
as functions of the temperature for spin S* = 1 films with thicknesses ranging from 
A'" = 1 to = 19 layers calculated with the exchange anisotropy using the same 
parameters as in figure 4. The reorientation temperatures can be read from 
the curve in the N-T plane, which approaches the bulk value with increasing layer 
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thickness. Similar results were obtained in reference jU] when using the single-ion 
anisotropy. To have a direct comparison of the reorientation temperatures calculated 
with the exchange anisotropy {D = 0.7) and the single-ion anisotropy {K2 = 1), we 
display figure 6 . Because the exchange anisotropy parameter was fitted to give the 
same result for the Curie temperature as a calculation with the anisotropy parameter 
for the monolayer, the corresponding reorientation temperatures practically coincide 
in this case, whereas the reorientation temperatures turn out to be slightly higher for 
the exchange anisotropy calculations with increasing film thickness. The saturation 
towards the bulk limit follows the same trend in both cases. 

reorientation temperature versus film tliicl<ness 




I > I > I > I > I > I > I > I > I > L- 

1 3 5 7 9 11 13 15 17 19 

number of layers 



Figure 6: Reorientation temperature as a function of the film thickness displayed 
for results with the exchange anisotropy and the Anderson- Callen treatment of the 
single- ion anisotropy (see reference |S] for the latter). 
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Reorientation angle for different film thicknesses 




temperature 



Figure 7: Average equilibrium reorientation angle for different film thicknesses 
(number of layers N=l,..., 5, 7, 9, 11, 15, 19) shown as a function of the temperature. 
This is a result of calculations using the exchange anisotropy and the parameters 
are the same as in the previous figures. 



In figure 7 we show the average reorientation angle ^o(-^, T) as a function of the 
temperature for films with increasing film thickness, where we define 

J_y-Af /qz\ 

^o(A^,r) = arctanf^^^. (58) 

The curves show the same saturation behaviour in the bulk limit as already seen in 
the previous figures. 
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4. Conclusions 

In the present paper, we have included the exchange anisotropy in our many-body 
Green's function description of thin ferromagnetic films. A few non-trivial changes 
in the general formalism had to be implemented because of additional momentum 
dependencies when the exchange anisotropy is included. The momentum depen- 
dence of the projector onto the zero eigenvalue space has been eliminated by an 
appropriate transformation. A novel feature is the appearance of complex eigenval- 
ues and eigenvectors of the non-symmetric matrix governing the equation of motion 
for the Green's functions when including the exchange anisotropy and the dipole 
coupling together. The complex eigenvalues occur quite naturally in the theory and 
have to be taken seriously. They are necessary to obtain correct results above the 
temperature where the magnetization (S^) drops below a certain value. 

The physical quantities calculated with the exchange anisotropy or with the 
single-ion anisotropy are fairly similar when the exchange anisotropy strength is 
fitted in such a way that it gives the same Curie temperature for a spin S — 1 
Heisenberg monolayer as that calculated with the single-ion anisotropy parameter 
of a certain strength, the rest of the parameters such as the exchange interaction 
strength and the dipole coupling strength being the same. 

In the present paper, all intra- and interlayer coupling parameters have been 
taken to be the same but the computer program is written in such a way that they 
can easily be chosen differently. We have also shown only examples for a monolayer 
with spins in the range between S — 1 and S — 6 and films with N layers for spin 
S — 1. It is only a question of computer time to make calculations for films with 
other {S, N) combinations. 
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Appendix A: Treating S > 1 

In this Appendix we sliow liow tlie regularity conditions can be deduced for general 
spin quantum numbers S and for multilayers. From the definitions ()12|) we see from 
dnSI) that 

TJX TDX 

The regularity conditions ()34|) can therefore be written for general m, n in the form 
for each layer i 

- 2ZiA'_lT; = At\TBt + AZ'U'B^ . (60) 

For the calculation of the correlations for higher spin we use equation (j25|) generalized 
to the multilayer case. We leave out the layer index i in all formulas which follow. 

- {{STiS'TS-) - -{^-fTi - <=-i) = 
l{A-r.1 + E ^ coth(^). (61) 

We express all correlation functions occuring in this equation in a standard form 
where all powers of are written to the left of the powers of 5*": 

C{m,n) = {{ST'{S~)''). (62) 

Then, with the relations (5-)"]_ = -n{S~)'' and 5" 5+ = S{S + 1) - -{S')\ 
we find that 

{{ST{S-)''S') = nC{m,n) + C{m+l,n) , 

{{ST{S-)''S+) = [S{S + l)-n{n-l))C{m,n-l)-{2n-l)C{m + l,n-l) 

-C(m + 2,n- 1) , 

{{STiS-rS-} = C{m,n + 1). (63) 

The commutators can also be expressed in terms of the C{m,n) using the binomial 
series 

A^Y — —nC (m, n) , 

A+f " = ([((5" - I)'" - {S')"')S-S+ + 2S'{S' - 1)"* + {n- l){n + 2S'){ST'] (5^)""^) 

m / ^ \ 

= 5(5 + l)E \ {-iyC{m-i,n-l) + {2n + m)C{m+l,n-l) 



m+l 



m + 1 



+ E ] {-iy^^C{m + 2-i,n-l) + n{n-l)C{m,n-l) 



i=2 



Ajr = ([(5^ + l)'"-(5^)'"](5-)"+^) = E ]C(m-z,n + l). (64) 
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Now by putting equation (j64j) into equation (j6Uj) the regularity conditions for all m 
and n can be written in terms of correlations defined in the standard form: 

r ™ / m \ 

2ZnC{m,n) = B''[S{S + j (-l)*C(m - n - 1) 

m+i [ rn+l\ 

+ (2n + m)C(m + l,n-l)+ ^ ( j (-l)*+^C(m + 2 - ^, n - 1) 

+n(n- l)C(m,n- 1) +5^^ \C{m-i,n + l) . (65) 



i=l 



For a given spin S*, this set of linear equations for the correlations has to be solved 
for all m + < 25* + 1. The solutions have to be put via equations (j63p together 
with (j64|) into equations (j6ip . thus leading to a set of 25* equations for the moments 
{{S^y) (p=l,. • • ,2S), which have to be solved self-consistently. The highest moment 
{{S^y^^'^) has been eliminated in favour of the lower ones through the relation 
UMsiS' - Ms) = 0. 
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